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1. Introduction 

The idea of volume independence in large N gauge theories goes back to Witten's original 
argument [1] that the path integral of N = °° Yang-Mills theory (YM) should be localized on a 
particular translational invariant field configuration, the master field. The translational invariance 
of the master field would imply the translational invariance of physical observables, 

{0) L= J=~ {@) L=0 . (1.1) 

Consequently, spacetime degrees of freedom become spurious in the description of large N physics. 

In particular, a matrix model would exist that describes the physics of pure SU (N) YM on M 4 , 
in the N — > oo limit. The possibility of such a formulation, originally proposed in [2], would allow 
a better analytical control over YM's complicated dynamics. On the numerical side, the absence 
of spacetime degrees of freedom would also allow faster simulations of the theory at larger N. 
This idea of large ,/V volume reduction can be applied to other gauge theories with a similar free 
parameter ,/V, or to the volume reduction of only a few spacetime directions (partial reduction). 

However, the reduction of spacetime degrees of freedom is not freely granted. In order for 
a SU (N) gauge theory to be independent of the volume of a particular (compact) direction, the 
center symmetry associated with it must be intact. Technically, the expectation values of the 
non-contractible holonomies £l x wrapping the reduced directions (a.k.a. reduced Polyakov loops) 
must be zero [2], 

(TrO x ) fix = J %(TrO x )=0, z € Z N . (1.2) 
Consider four-dimensional SU (N) lattice gauge theory with the standard Wilson action, 

Sw(t/) = -2NX-' £ £ ReTr(U^ x U v x+fl Ul x ^Ul x ), (1.3) 

where X = Ng 2 is the lattice 't Hooft coupling. In this theory, large ,/V volume independence only 
occurs when the volume of a reduced direction is larger than a critical value, i.e. L > L c , where a 
Z^-symmetric confining phase appears [3]. For L < L c , the Z;v symmetry breaks and the expecta- 
tion values of reduced Polyakov loops become non-zero (deconfined phase). Hence large ,/V volume 
reduction does not hold in general. In order for the volume reduction to hold for arbitrarily small 
volumes, the lattice action (1.3) needs to be modified in a way that stabilizes the Z/v symmetry. 

Consider pure YM compactified on R 3 xS\ and let L be the volume of the compact direction 
(Fig.l). At small L, the symmetry is broken along the compact direction. The reason behind it 
resides in the fact that the effective potential of reduced Polyakov loops [4] , 

i W2\ j 

MA) = - j3 I I ^4> Tr ^l 2 (L4) 

is minimized at maximal traces: (Tr £l n ) ^ 0, Vn. The Z^ symmetry would be restored if the sign 
of (1.4) was flipped, a situation in which V e s would minimize at vanishing traces. Such a stabiliza- 
tion of the Zn symmetry can be achieved by either adding massive adjoint fermions with periodic 
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Figure 1: Reduced Polyakov loop wrapping the compact direction of M? x S . 

boundary conditions along the reduced directions [5], or by adding double-trace deformations that 
counteract the Zjv-breaking character of the effective potential [6], 

l [N/2] 

5 d YM(f/)=5 w (f/) + 7 3 £ £ a n \TrD. n x \ 2 , (1.5) 

xeR 3 n=\ 

known as deformed YM (dYM). For sufficiently large deformation parameters a n the Zn symmetry 
is restored for all L, making dYM theory on R 3 x S 1 equivalent to pure YM on M 4 in the N — > °° 
limit (where deformations become irrelevant). See also M. Unsal's talk in these Proceeedings [7]. 

In this talk we summarize the results of our non-perturbative study [8] of the phase structure of 
dYM on R 3 x S 1 (and related question of large N volume independence), for which we performed 
numerical simulations on a lattice with a fully reduced direction (L = 1). 



2. Monte Carlo algorithm 

The numerical simulation of SU (N) lattice gauge theories with the standard Wilson action 
(1.3) can be performed efficiently using the Cabibbo-Marinari pseudo-heatbath algorithm [9]. A 
Metropolis algorithm could also be used, but it is typically slower, less ergodic, needs tuning for op- 
timal acceptance rates, and has larger autocorrelation times. Pseudo-heatbath algorithms are faster 
and require no tuning, but they can only be applied to lattice actions that are linear with respect to 
each link variable. If that condition is satisfied, the probability distribution for a link U reduces to 

p(U)=e R ^ v ' tu \ (2.1) 

where V is the relevant contribution from neighboring links, which must not include U explicitly. 

In the case of dYM, a pseudo-heatbath algorithm cannot be applied directly, since the deforma- 
tion terms in (1.5) are not linear in the link variables. However, by introducing a sufficient number 
of auxiliary free scalar lattice fields M,- in the dYM partition function, and performing adequate 
Hubbard-Stratonovich transformations, M,- 1-> M,(M), 

Z dYM = JDU e- s *™(u)^ j d M. d Mle- l 2 T ^M) (2 .2) 



M^um J DU e - hMM ^M) e - S ^U) e K^ { V { UMM^O.)^ (2J) 

we are able to linearize the dYM action (1.5) with respect to the link variables along the compact 
direction, U = Q. X . Consequently, we can construct a useful pseudo-heatbath algorithm for lattice 
dYM [10]. This trick was inspired on a pseudo-heatbath algorithm constructed in a similar way by 
Fabricius and Haan for the twisted Eguchi-Kawai reduced models [11]. 
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The auxiliary fields M, = (R n ,xjQn$) and respective Hubbard-Stratonovich transformations 
that linearize (1.5) are given by 1 



R 



7/.X 



V2NL-^a n (R, hX - - tr 1)) , 1 < n < [N /2\ , 

tT(R n ,x% 

OH = V2NL-ia n (eg - (efc _1) nj + Q^ m )), l<m<n< \Nj2\ . 
(jtl = V2NL-^a n (eS - tr(/?„ iX )a+), 2 < n < m = [N /2\ . 



X(m) 



(2.4) 
(2.5) 
(2.6) 
(2.7) 



For the update of the links along the compact direction, £2 X , the contribution V = V X coming 
from the neighboring "staples" and deformations is given by 

3 



V x = 2NX' 



I 

v=l 



L { J+L^)+2NL-'f x , 



(2.8) 



where £ x J and are the usual backward and forward "staple" contributions at x, and f x is the 
contribution coming from the deformations, encoded as a function of the auxiliary fields, 



[N/2\ ( n-\ 

J x — a\Q\. x + 2- a n\\2n,y. + Un,x Qn,x + ^{R n ,x)Un,x 

n=-1 \ m=l , 



(2.9) 



We tested the pseudo-heatbath algorithm by using it in simulations of SU (5) dYM on a 10 3 1 
lattice, for several couplings A -1 and deformation parameters (01,02). which we compared with 
simulations of the same parameters using a Cabibbo-Marinari-Metropolis algorithm [12]. Expecta- 
tion values of the plaquettes coincide in both algorithms for all cases, which confirms the validity 
of the pseudo-heatbath algorithm. In the pseudo-heatbath case, however, there was a significant 
improvement over Metropolis in terms of equilibration and autocorrelation times (see Fig. 2). 
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Figure 2: Estimator of the autocorrelation function (C) vs. CPU time (fcpu) i n simulations of SU (5) dYM 
theory on a 10 3 1 lattice, with A -1 =0.5 and (01,02) = (0.20,0.05) (left), and the corresponding estimator of 
the integrated autocorrelation time (Tj nt ) vs. Monte Carlo time (t) (right). We compared a Cabibbo-Marinari- 
Metropolis algorithm for dYM (red) with the pseudo-heatbath algorithm described in the text (blue). 



'To we denote the normalized trace by tr = ATr. 
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3. Phase structure 

Using the pseudo-heatbath algorithm described in the previous section, we simulated SU (N) 
dYM on a 6 3 1 lattice, for N = 4, 5. 2 We determined the phase diagrams of these theories (Figs.3- 
4) as functions of the coupling A -1 and of the single- winding deformation parameter a\. For that 
purpose we calculated the expectation value of reduced Polyakov loops, (tr Q£), which are order 
parameters of the Zjy symmetry associated with the reduced direction, and its subgroups. 

For all cases, at small a\, there is a bulk transition between a Z^-symmetric strongly-coupled 
bulk phase (left) and a Z^-broken weakly-coupled deconfined phase (right). The hysteresis region 
(shaded area) suggests that this is a first-order transition. The breaking of the Zn symmetry in this 
regime of the deformation parameters indicates, as expected, that the contribution of the deforma- 
tion terms in the dYM action is negligible and is not sufficient to restore the symmetry at small 
volumes. 

At large fixed a\ there is a smoother bulk transition that occurs at a fixed value of the lattice 
coupling, At A -1 < X~ , the phase is the same Z^-symmetric strongly-coupled bulk phase. 
At A. -1 > A c , however, the weakly-coupled phase is not always deconfming. If we consider a 
large value of the double- winding deformation parameter, e.g. t?2 = 0.1, the weakly-coupled phase 
is Z;v-symmetric, i.e. confining. This indicates that the double-trace deformations are large enough 
to counteract the Z^-breaking character of (1.4), thus preserving the Z^ symmetry at all couplings. 

At large a\ but vanishing ai, however, the Zm symmetry is partially broken to a subgroup 
(Z2 for N = 4, and Zi for N = 5). This is due to the fact that ai is not large enough to preserve 
the subgroup Z p C Z^ of which (tr Q. 2 } is an order parameter. Therefore, the transition Z^ — > Z p 
occurs. 

In sum, the Z^ symmetry in dYM is preserved at all couplings only when each of the deforma- 
tion parameters a n is sufficiently large. In this situation, a large N orbifold equivalence should hold, 
and dYM on M 3 x 5 1 should reproduce ordinary YM on M 4 up to 0(1 /N 2 ) corrections. With this 
in mind, we also simulated SU (N) YM with the standard Wilson action on a 6 4 lattice. We located 
the bulk transition for ,/V = 4,5 and plotted it as a vertical magenta line in Figs. 3 -4. We observe 
that the bulk transitions of dYM on a 6 3 1 lattice and of YM on a 6 4 lattice coincide when the Z^ 
symmetry is intact in the weakly-coupled phase; but they do not coincide when the Z^ symmetry is 
(partially) broken there. This result indicates that the expected large ,/V equivalence between dYM 
on M 3 x S l and YM on M 4 may hold quite accurately even for small N. 

There should also be a phase transition between the regime of small deformations, a\ < a\ >c 
(where the Zjv symmetry is broken), and the regime of large deformations, a\ > a\ yC (where Zn 
symmetry is intact and volume reduction holds). Unsal and Yaffe suggest that the critical values of 
the deformation parameters are given by a, hC = 4/n 2 n 2 [6]. Such a choice corresponds to adding 
— 2V e ff to the standard Wilson action (1.3). Perturbatively, the effective potential of dYM would be 
—V e s, which minimizes at vanishing traces of the the reduced Polyakov loops, (Tr £1£) = 0. This 
prediction is represented in Figs. 3-4 by an horizontal green line. 

In our simulations we observe a large hysteresis associated with the transition between small 
and large values of a\, which suggests that it must be strongly first-order. The transition lines at 

2 For these gauge groups, dYM requires a pair of independent double-trace deformations, given by single- and 
double-winding reduced Polyakov loops, and weighted respectively by the free parameters a\ and a*i- 
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Figure 3: Phase diagrams on the (Wi) plane of 51/(4) dYM on a 6 1 lattice, for vanishing double- 
winding deformation (left) and large double-winding deformation (right). 
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Figure 4: Phase diagrams on the (A l ,a\) plane of 51/(5) dYM on a 6 3 1 lattice, for vanishing double- 
winding deformation (left) and large double-winding deformation (right). 



asymptotic values of the lattice coupling roughly approach the value ax >c ~ 4/n 2 , which supports 
the prediction of Unsal and Yaffe. 



4. Conclusions 

We constructed an efficient Monte Carlo algorithm for SU (N) dYM on R 3 x S 1 that allows the 
direct use of the Cabibbo-Marinari algorithm for updates of the link variables. This algorithm has 
better equilibration and autocorrelation times than an optimized Metropolis algorithm for dYM. 
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We simulated SU(N) dYM on a 6 3 1 lattice for /V = 4,5 and mapped their phase diagrams 
with respect to the expectation values of (multi-winding) reduced Polyakov loops, which are order 
parameters of the Z# symmetry associated with the fully reduced direction (and its subgroups). 

Lattice dYM possesses a rich phase structure, with phases similar to the ones discussed in 
[13]. These include confining, deconfining and partially confining phases. 

At vanishing values of the double- winding deformation parameter, ci2, the center symmetry 
breaks partially to a subgroup of Zn, because (tr £2^) acquires a non-zero value. Only for sufficiently 
large values of both a\ and ci2 is the symmetry fully preserved at all couplings. 

When the Zn symmetry is intact at all couplings, lattice dYM on M 3 x S l seems to reproduce 
ordinary lattice YM on M 4 . This is suggested by the fact that the critical couplings of the bulk 
transition in both theories coincide quite accurately. 
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